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ABSTRACT 

We investigate dynamical response on size and velocity dispersion to mass loss by super- 
novae in formation of two-component spherical galaxies composed of baryon and dark matter. 
Three-dimensional deprojected de Vaucouleurs-like and exponential-like profiles for baryon, 
embedded in truncated singular isothermal and homogeneous profiles for dark matter, are con- 
sidered. As a more realistic case, we also consider a dark matter profile proposed by Navarro, 
Frenk & White. For simplicity we assume that dark matter distribution is not affected by mass 
loss and that the change of baryonic matter distribution is homologous. We found that the 
degree of the response depends on the fraction of dark matter in the region where baryon is 
distributed, so that dwarf spheroidal galaxies would be affected even in a dark halo if they 
are formed by galaxy mergers in the envelope of the dark halo. Our results suggest that this 
scenario, combined with dynamical response, would make not only the observed trends but 
the dispersed characteristics of dwarf spheroidals. 

Key words: galaxies: dwarf - galaxies: elliptical and lenticular, cD - galaxies: evolution - 
galaxies: formation - galaxies: haloes - large-scale structure of the universe 



1 INTRODUCTION 



Elliptical galaxies have provided important information on the formation and evolution of galaxies. So far some scaling relations among 
photometric properties and global struc tural parameters have been observed, such as the colour-magnitude relatio n (e.g ., Baum 195S| ), the 



& Watanabe 



velocity disp ersion -magnitude relation ( Faber & Jackson 1976 ), and the surface brightness-size relation ( Kormendy 



1977 



Kodaira, Okamura 



1983 ). These relations are regarded as clues to our understanding of galaxy formation. 



It is widely believed that heating up and sweeping out the galactic gas by multiple supernova (SN) explosions play a key role on 
the formation of dwarf elliptical/spheroidal galaxies. In a traditional framework of monolithic cloud collaps e scenario fo r formation of 
ellipticals, s uch multiple SNe cause a galac ti c wind to halt subsequent s tar formation and chemical enrichment ( Larson 1969 ; l^keuchi 



Saito 



1979| ; [Arimoto & Yoshii||l986 , 



1987; Kodama & Arimoto 



Hills 



1977 



1997 ). Since a large amount of gas is expelled by the galactic wind, the 



198C; Mathieu 



1983 



Vadei 



1980. Yoshii & Arimoto (19871 found that 



self-gravitating system expands after blowing up the wind (e.g., 
many observed properties of ellipticals can be reproduced by the galactic wind model with the evolutionary population synthesis technique. 
They noticed that theoretical metallicity of galaxies, t o be compared with observation, should be an average of metallicities of constituent 
stars weighted by their luminosities. On the other hand, Dekel & Silk (19861 considered a model in which galaxies are embedded in dominant 
dark haloes. They claimed that structural properties of luminous matter would not be affected by the gas removal if dark matter dominates 
the gravitational potential. Their model can also reproduce the observed relations, although they did not take into account the evolutionary 
population synthesis. Recently, Zhao (2002) argued the dynamical limit on occurrence of galactic wind given by the existence of sufficiently 
dense dark halo and the absence of intergalactic/escaping globular clusters, and concluded that only about 4 per cent of the total mass of the 
galaxy should be lost by the wind. 

In contrast to the monolithic cloud collapse scenario, recent studies on the cosmological structure formation suggest the hierarchical 
clustering scenario based on a cold dark matter (CDM) model. In this scenario, the CDM dominates the Universe gravitationally and objects 
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or dark haloes continuously cluster and merge together, so that larger objects form via mergers of smaller objects. Thus, simple galactic 
wind models must be modified in order to be consistent with the cosmological structure formation. Based on the hierarchical clustering 
scenario, semi-analytic models (SAMs) of galaxy formation have been developed (e.g., Kauffmann, White & Guiderdoni 1993; |Cole et aF 



1994, 200C 



Somerville & Primack 



1999; Nagashima, Gouda & Sugiurt 1999; Nagashima et al 



2001 



2002). These models well reproduce 



many statistical observables of galaxies such as luminosity function, colour distribution, and gas fraction. 



Kauffmann & Chariot 



( |1998[ ) and 

Nagashima & Gouda (2001 ) have shown that their SAMs taking into account metal enrichment can reproduce the observed colour-magnitude 



relation of elliptical galaxies in clusters of galaxies and that this relation reflects the mass-metallicity relation, as pointed out by using the 
galactic wind model (Kodama & Arimoto 1997). Thus, as a next step, it is valuable to investigat e not only the photometric properties but 
structural properties such as size and velocity dispersion. Moreover, Okamoto & Nagashima ( 2001a ^ found that radial gradients of fractional 
number, colour, and star formation rate of ellipticals in clusters are also reproduced by using the SAM associated with an TV-body simulation. 

Because a large amount of mass loss might occur during the formation of ellipticals, in both scenarios of monolithic cloud collapse and 
hierarchical clustering, it is important to construct a formalism on the dynamical response in two-component galaxies consisting of baryon 
and dark matter. While self-consistent equilibrium models of such two-component elliptical galaxies have been analysed by Yoshii & Saio 
( 1987) and Ciotti & Pellegrini ( 1992), we focus on the dynamical response of baryon within a dark matter halo. In this paper we consider 
simple models of two-component galaxies in which the baryon having either de Vaucouleurs or exponential profile is embedded in the dark 
matter having either isothermal or homogeneous profile. Moreover, we also consider the so-called NFW profile of dark matter that is recently 
suggested by Af-body simulations (Navarro, Frenk & White 1997). 

This paper is outlined as follows. In section 2, general formulation for the dynamical response on size and velocity dispersion to mass 
loss is presented. In section 3 the dynamical response for particular choices of density distributions mentioned above is investigated. Section 
4 summarizes the results of this paper. The details of models and derivations are given in Appendices. 



2 EQUATIONS OF DYNAMICAL RESPONSE IN TWO-COMPONENT GALAXIES 



We assume that the velocity fields of gas, stars and dark matter are non-rotating, homogeneous and isotropic characterized by a single 
parameter of velocity dispersion a, and that the density distribution of baryon changes homologously. For simplicity, it is also assumed that 
the density distribution of dark matter does not change during the mass loss. This assumption will be justified for global structural parameters 
of dark matter because the total mass of dark matter dominates over baryon by a factor of 10 or so. Although the following results will be 
slightly modified in a dense region of baryonic matter because of non-negligible dynamical response on the local dark matter distribution, the 
resulting change of baryonic matter should be intermediate between the baryon-dominated and dark matter-dominated cases. We will study 
the dynamical response in such limiting cases in section 2.2. 

The response depends on the timescale of gas removal between two extreme cases of adiabatic (slow) and instantaneous (rapid) gas 
removal. When the timescale is longer than the dynamical timescale of the system, it changes adiabatically in quasi-equilibrium. In this 
case the total change is a sum of consecutive infinitesimal changes. On the other hand, when the timescale of gas removal is much shorter 
than the dynamical timescale, the system drastically changes and it is possible that our assumption of homologous change might be broken. 
Under realistic circumstances, as hydrodynamical simulations suggest, these timescales are almost similar (e.g., Mori, Yoshii, Tsujimoto & 
Nomoto 



1997; Mori, Yoshii & Nomoto 



1999; Mori, Ferrara & Madai 2002). Therefore, keeping in mind that the adiabatic gas removal is 



more consistent with the situation under our assumptions, we consider both adiabatic and instantaneous gas removal below. 



2.1 Formulation 

A kinetic energy T due to velocity dispersion and a gravitational potential energy W of baryonic component in a dark matter halo are written 
as 

T = ^PPtrla^ (1) 

W = -aplr^ - bpbPdrtraf{z), (2) 

respectively, where p, a and b are constants dependent on density distribution, pt and pd are characteristic densities of baryon and dark 
matter, respectively, a is velocity dispersion of baryon, rt and rd are characteristic radii of baryon and dark matter, respectively, z is a radius 
rt normalized by r^, z = rt/rd, and f{z) is a function dependent on density distribution. Note that variables related to dark matter, with 
subscript d, do not change during the mass loss. From the virial theorem, we obtain 

pa^ = apir^ + hpdrif{zi), (3) 

paj = apfrj + bpdrjf{zf), (4) 

where subscripts i and / stand for the initial and final states, that is, before and after the mass loss, respectively, for the baryonic component. 
Just after the mass loss, only the baryon density changes from pi to p'^ without any changes of other variables. Then, conserving its energy 
and mass, the system virializes again quickly to the final state. 
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-ppiTiCJi - ap, ri - bp,pdrira,f{zt) = -ppfTfOf - ap/ry - bpfpdrfra.T(zf). (5) 
After the mass loss the total mass of baryon also conserves, then p'^rf — Pfr'j. By using this relation, the above equation is transformed to 



^ a^) - apfvj (l - ^) - bp^rllfizf) - f{z,)] = 0. 
Using the virial relation, ct^ and cr/ can be eliminated as 

^^'r[/(^/)-/(^.)] = 0. 



l + EL 



(!)■ 



ap.zf 



When the gas removal occurs instantaneously, the density at the final state is given by the above equation, 

Vf _ 1 + b[f{zf) - f{zi)]/ayaf 



(6) 



(7) 



(8) 



Vi (zf/zi)^{2zf/zi-l) 

where yi,f = pij/pd, respectively. On the other hand, when the gas removal is sufficiently slow, the total change is a sum of consecutive 
infinitesimal changes, for which y f = yi + dy and Zf — Zi + dz. Thus, linearizing the above equation, we obtain 



dy 
dz 



,V , b df{z) 



(9) 



z az^ dz ' 

where subscript i is omitted for simplicity. This corresponds to the case of adiabatic change. The first term in the above equation is obtained 
for the self-gravitational potential of baryon. The second term gives the modification for baryon embedded in dark matter. Solving this 
differential equation, we obtain 



C 

y = — + <i{z), 

where 



t^'J^dt, 
dt 



(10) 



(11) 



and C is an integration constant. Since C is an adiabatic constant, the value does not change during the mass loss. Thus we obtain the density 
at the final state. 



and 



C/z*j+q{z!) ^ l_ 



{yi - <l{zi))-^Jrq{zf ) 



C ^ z^[y- q{z)] = zf[yi - q{zi)]. 

From equations (pi) and (W), the velocity dispersion at the final state is obtained as 



5L 



yfz} + bf{zf)/a 



1 1/2 



yiZi +bfiz^)/a 



(12) 



(13) 



(14) 



In the next section, we derive the dependence of y on z by specifying the function f{z) in each case of density distributions of baryon 
and dark matter. Before studying such specific cases, we see the limiting cases of no dark matter and negligible baryonic matter in the next 
subsection. 



2.2 Limiting cases 

First, w e consider the case of self -gravitating system with no dark halo, that is, pd/pi 
result in 



0. The final form is expected to be the same as the 



Yoshii & Arimotc (1987). For instantaneous gas removal, we obtain from equation i 



Mf/M, 



1' 



(15) 



n 2Mf/Mi 

where Mj /Mi = pfr'^f / pir^- This is the well-known result that the system becomes unbound if a half of the mass is removed. On the other 
hand, for adiabatic gas removal, the second term in equation vanishes and then 

^4 



El 

pi 



The velocity dispersion is also obtained from equation ([14 

1/2 



ZL 



(16) 



(17) 
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Next, we consider the opposite case, pi/ pd — > 0. As pointed out by Dekel & Silk (1986), the size and velocity dispersion of galaxies 
do not change during the mass loss because the gravitational potential is dominated by dark matter. Since the third term in equation 
overwhelms other terms, z does not change against the change of y. Thus the velocity dispersion also does not change as obtained by 
equation (p^. This case corresponds to that of §III-c in Dekel & Silk ( 198(: ). 

Therefore, we confirm that our formulation is applicable to the limiting cases of no dark matter and negligible baryonic matter. 



3 DYNAMICAL RESPONSE ON SIZE AND VELOCITY DISPERSION IN SPECIFIC CASES 



In this section, we consider some specific density distributions. As for the baryonic component, we adopt the following two distributions 
One is the Jaffe model (laffe 1983) that approximates the de Vaucouleurs profile when projected, 

'^Pbrt 



P{r 



(18) 



where pb and rj, are the characteristic density and radius, respectively. Another is a model that approximates the exponential profile when 
projected. 



p{r) 



Pb\ 



— exp 

r 



(19) 



The details of comparisons of these approximate profiles with the de Vaucouleurs and exponential profiles are given in Appendix A. As for 
the dark matter halo, we consider the following two distributions for truncated isothermal and homogeneous spheres. 



p{r) = p^^e{rd^r) 



p{r) 



pdO(rd 



(isothermal) 
(homogeneous) 



(20) 
(21) 



respectively, where pd is a characteristic density of dark matter and 9{x) is the Heaviside step function. Finally we also consider the NFW 
profile (iNavarro, Frenk & White||l997|), 



p(r) = pdC 



cr 








Td 





(22) 



where c is the concentration parameter. Although ^^ukushige & Makino (1997, 2001 ) and Moore et al. (1998) suggest a modified NFW profile 
with steeper slope in the inner core, p oc r^^'^, we consider only the original NFW profile. 

In a hierarchical clustering scenario, the isothermal profile should be applied to mergers and starbursts occurring near the centre of 
dark halo. This situation is realized when the mergers with satellite galaxies are caused by the dynamical friction. On the other hand, when 
satellite-satellite mergers occur in the envelope of dark halo by random collisions, the homogeneous profile should be applied. Recent high 



resolution A*'-body simulations suggest that many subhaloes survive and are not disrupted (e.g., Okamoto & Habe 



19% 



Moore et al. 



1999 



Ghigna et al. 2000). If their internal structures are also maintained, the isothermal profile should be applied even in the case of mergers in the 
envelope of dark halo. Moreover, outer parts of individual subhaloes near the centre of their host halo would be stripped by the tidal force 
from the host halo, then the mass ratio of baryon to dark matter would be high. We also consider such a case below. 



3.1 de Vaucouleurs-like distribution of baryon in an isothermal dark halo 

In this subsection, we investigate the dynamical response to mass loss in the case of isothermal dark halo. As already mentioned, we adopt 
the Jaffe model which well reproduces the deprojected three-dimensional de Vaucouleurs profile. The corresponding mass profile for baryon 
is 



Mi,(r) = 



IGTrpbr^r 



rd 



r + rb 

and for dark matter, 

Md{r) = ^TTpdT-d — 

ITd \rd 

In the following, we use a parameter m that stands for the mass ratio of baryon to dark matter, 

Mb(r oo) 
Md{r <x)' 

From the self- and interaction-potential energies of baryon, we obtain the coefficients a and b in equation (^ as follows, 

a = 1287r^G, and b = 64-k^G, 
and the function f{z) is 



(23) 



(24) 



(25) 



(26) 
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ln(l + ^) 



In 



1 + i 

z 



(27) 



Substituting these quantities into equation ( pi] ) for adiabatic gas removal and into equation (^) for instantaneous gas removal, we obtain the 
dynamical response on size to gas removal. Moreover, using equation (jl^), we also obtain the response on velocity dispersion. Because the 
corresponding equations and solutions are complicated, we summarize those in Appendix ^ 

We examine how the dynamical response depends on the ratio of initial mass of baryon relative to dark matter halo nii = Mi /Md and 
the ratio of initial radius of baryonic component relative to dark matter halo Zi = Vi/rd- In the following we consider two representative 
values of nii = 0.2 and 0.05, and Zi — 0.2 and 0.05, respectively. In a standard cosmological model, the baryon fraction is about 0.1, so 
rrii ~ 0.1 is expected if gas cooling is sufficiently effective. Thus the two values of correspond to two extreme cases in realistic situations. 
The size of disk galaxies is usually considered to be a radius at which the disks are supported by their rotation under an assumption that their 
specific angular momenta conserve while the gas shrinks due to the cooling. The initial value of the specific angular momentum before the 
cooling is given by the so-called dimensionless spin parameter A, which distributes log-normally around 0.05, derived by the tidal torque due 



1984; Catelan & Theuns 



1996a ,:; Nagashima & Gouda 1997). Because an 



to the large-scale density inhomogeneity in the Universe (White 1 
effective radius of disk is ~ Arvir (Fall 1983) and the size of elliptical galaxies is not significantly different from the disk size, the values 
of Zi = 0.05 and 0.2 are realistic. If outer parts of individual subhaloes are tidally stripped by their host halo, even larger values might be 
realized. 

In the left panels of Figure |l|, we plot the responses on size (upper panel) and velocity dispersion (lower panel) for rrii = 0.2. The 
thick solid and dashed lines indicate the results of Zi — 0.2 and 0.05, respectively, for adiabatic gas removal. In the latter case of Zi = 0.05, 
because the baryonic component is compact and gravitationally dominates in the central region, the changes of size and velocity dispersion 
are larger than in the former case of = 0.2. In the right panels, we show the results of rrii = 0.05. The changes are smaller when compared 
to those for nii = 0.2 because of the small mass of baryonic component. Note that the changes for {rrii, Zi) — (0.2, 0.2) and (0.05, 0.05) are 
very similar to each other. The reason is, as pointed out by Navarro, Eke & Frenk (I99t ), that the changes can be measured by the value of 
rrii/zi. If rrii/zi is sufficiently small as shown by the thick solid line in the right panels (nii/zi — 0.25), the changes are negligible because 
the mass of dark matter dominates over the baryonic component. In such a case the assumption of dominant halo by Dekel & Silk (1986) can 
be justified. 

For reference, we show the results for instantaneous gas removal by thin lines. As well known, the changes are greater than those for 
adiabatic gas removal. This tendency holds even in the dark matter-dominated case. Note that our results for this case are obtained under the 
assumption of homologous change, which might not be justified because some masses become unbound and leave the system. Actually, by 
using A'^-body simulations, Navarro, Eke & Frenk (1996) showed that some particles initially in the central region of a dark halo become 
unbound when a part of mass is instantaneously removed. 



3.2 Dependence on dark matter distribution 



In this subsection, we examine how the dynamical response depends on the density distribution of dark matter component, while the density 
distribution of baryonic component is fixed to the Jaffe model. 

First we consider that the dark matter distributes homogeneously. The mass of dark matter within radius r is 



Md(r) 



An 3 



^{r - rd) 



Then the coefficient b is 

. 64 2 
h=-. G, 

and the function f[z) is 



1-22 ^ ln(l + z) 



z^lnil 



In the same manner as before, we obtain the dynamical response on size and velocity dispersion for this case. 
Next we consider the NEW profile. The mass within radius r is 



Mdir) = Aixpdri 



In 



V Td 



fjrd 



1 -I- crjrd 



The coefficient 6 is 

6 = 647r^(7, 

and the function /(z) is 

j(^) ^ £M££l + 1 [in(l - cz) In(cz) + Li2(c2)] . 
\ — cz z 



(28) 



(29) 



(30) 



(31) 



(32) 



(33) 
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Figure 1. Dynamical response of the baiyonic component of deprojected de Vaucouleurs profile embedded in an isothermal halo, for various values of initial 
mass mi = Mi/M^ and size Zi = ri/r^. The horizontal axis is the ratio of final to initial baryon masses. Upper and lower panels indicate the changes of 
size and velocity dispersion, respectively, for m,i = 0.2 (left panel) and 0.05 (right panel). The thick solid and thick dashed lines in each panel show the results 
of Zi = 0.2 and 0.05, respectively, for adiabatic gas removal. The thin solid and thin dashed lines show those results for instantaneous gas removal. Note that 
the scales of vertical axes of left and right panels are different. 



where hi2{z) is the dilogarithm. In defining the parameter m in equation (|25|), we truncate Md{r) at r = r^, because it logarithmically 
diverges at r cxj. Otherwise we do not assume any truncation in the analysis. 

The results for adiabatic gas removal are shown in Figure |^ The initial condition is m,; = 0.2 and Zi = 0.05 for all cases and c = 10 for 
the NFW profile. The solid, dashed and dot-dashed lines indicate the isothermal, homogeneous and NEW profiles, respectively. In the case 
of the homogeneous dark halo, we find that the changes of size and velocity dispersion are larger than those in the case of the isothermal halo 
and the influence from dark halo is weaker. This would reflects that the concentration of dark matter is lower than that of the isothermal halo 
and that the gravitational potential energy from dark matter is smaller, as indicated by the value of the coefficient b. This suggests that the 
size and velocity dispersion can be significantly affected by mass loss even if surrounded by dark matter, particularly when galaxy mergers 
frequently occur in the envelope of dark halo. 

In the case of the NFW halo with c = 10, the responses are very similar to those for the isothermal halo. This is because these haloes 
have much the same density distributions. 



3.3 Dependence on baryonic matter distribution 



In this subsection, we examine how the dynamical response depends on the density distribution of baryon within a given halo. We consider 
the exponential-like distribution in the isothermal halo. The mass within r is 



Mi{r) 



3 Tre-^ 



,5/2 



(2^ + 3) exp(-^)+3V^erf(^y^ 



(34) 



where we replace gexp by g for simplicity. Then the coefficients a and b are 
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Figure 2. Dynamical response of the deprojected de Vaucouleurs component within a dark halo with three different profiles. The solid, dashed and dot-dashed 
lines denote the isothermal, homogeneous and NFW profiles, respectively. Shown are the responses on size (left panel) and velocity dispersion (right panel) 
only for adiabatic gas removal for the initial condition of rrii = 0.2 and Zi = 0.05. 



a — — ^ — G, and b = — G, 



and the function f{z) is 
2g exp{-g/z) 



+ 



29 



23 



erf (yi) + In 64 + 3 In (^) } + S^iGig ( 



3/2,3/2 
1/2,1/2,1 I ' 



(35) 



(36) 



where 7 is the Euler constant and GJ^" is the Meijer G function (this symbol should not be confused with the gravitational constant G). 

The results for adiabatic gas removal are shown in Figure ^. Although we have calculated the exponential profile of baryonic component 
in the homogeneous halo (Appendix B4), we do not show this case in the figure, because the response is very similar to that for the de 
Vaucouleurs profile in the homogeneous halo. The dependence on the density distribution of baryonic component is weaker compared to that 
of dark matter component, which suggests that the response mainly depends on the density distribution of dark matter. 



4 SUMMARY 

We analyzed the dynamical response on structural parameters of elliptical galaxies to mass loss. While in the previous analyses only two 



limiting cases of pure self-gravitating baryonic system and dark matter-dominated system were considered (e.g., Yoshii & Arimoto 



1987 



Dekel & Silk 1 986 ), we extended the formalism to the intermediate case between the two. Three-dimensional density distributions of baryonic 
matter considered in this paper are the Jaffe model that approximates the de Vaucouleurs profile and a model that approximates the exponential 
profile in two dimensions. As for the density distribution of dark matter, a truncated isothermal sphere and a homogeneous sphere are 
considered. As a more realistic case, we also considered the NFW dark matter halo which surrounds the baryonic component having de 
Vaucouleurs profile. Thus we found, from five combinations of baryon and dark matter, that the response strongly depends on the central 
concentration of dark matter, reflecting the fraction of baryon relative to dark matter in the central region. In contrast to the simple assumption 



by Dekel & Silk (198( ) that structural parameters are not affected by mass loss, size and velocity dispersion can change even in dark matter, 
if baryon has sufficiently shrunk and become dense compared to dark matter at the onset of mass loss. 

In this paper we considered both adiabatic and instantaneous gas removal. The adiabatic removal is justified if the timescale of gas 
removal is sufficiently slow compared to the dynamical timescale of the system. On the other hand, in the case of instantaneous gas removal, 
the effects of gas removal are stronger. For example, when more than a half of the mass is removed instantaneously, the pure self-gravitating 
system becomes unbound. Thus, our results for adiabatic removal place a lower limit on the dynamical response to mass loss. In order to 
check the opposite case, we also derived the results for instantaneous gas removal. As expected, the changes of size and velocity dispersion 
are larger than those for adiabatic gas removal. Such drastic changes may not validate our simple assumptions of homologous expansion of 
baryon and static distribution of dark matter. 
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Figure 3. Dynamical response of the baryonic component within the isothermal dark halo. The solid and dashed lines denote the de Vaucouleurs and expo- 
nential profiles for the baryonic component, respectively. Shown are the responses on size (left panel) and velocity dispersion (right panel) only for adiabatic 
gas removal for the initial condition of rm = 0.2 and Zi = 0.05. 



So far some scaling relations such as the velocity dispersion-magnitude relation have been discussed considering the dynamical response 
to mass loss in the context of monolithic cloud collapse scenario for formation of elliptical galaxies, However, since the resuting structural 
change of baryonic component is found to depend significantly on the density distribution of dark matter, such scaling relations in the hier- 
archical clustering scenario should necessarily be reconsidered. Our results suggest that the dispersed properties of dwarf/compact ellipticals 
would be well reproduced by a variety of their formation histories and environment, when taking into account mass-dependent gas removal. 

In the hierarcliical clustering scenario, galaxy formation processes are rather complicated. There are many physical processes such as 
radiative gas cooling and galaxy merger as well as star formation and SN feedback. It should be noted that the physical mechanism of galaxy 
merger would be dynamical friction and/or random collision, while the mechanism of dark halo merger is growth of density fluctuation. In 
order to investigate the observed properties of galaxies, we need to fully incorporate the effects of dynamical response into realistic models 
of galaxy formation like our SAMs. We will discuss various observables of ellipticals in more realistic situations in a forthcoming paper. 
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APPENDIX A: SIMPLE MODELS OF DENSITY DISTRIBUTION 
Al Approximate density distribution for tlie de Vaucouleurs r^^* profile 

Most of elliptical galaxies have the so-called de Vaucouleurs r^'''* luminosity profile. 



fJ.{R) — Me exp 



(Al) 



where gdv ~ 7.67, R is the distance from the centre of luminosity distribution, and R^ is the effective half-light radius, at which the surface 
brightness reaches /le. In this paper, assuming that mass traces luminosity, we regard the surface brightness profile fi{R) as a surface mass 
density profile. 

Since the de Vaucouleurs r^^* profile is two-dimensional projected density distribution, we need to reconstruct three-dimensional 
density distribution. A direct procedure is to use the Abel inversion between a surface density profile and a three-dimensional density 
distribution p(r). 



KR) = 2 



p(r)r 



7?2 



I dr. 



(A2) 
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(a) Surface density profiles 
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(b) Surface mass profiles 
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Figure Al. (a) Surface density plotted against R. Ths solid, dashed, dot-dashed, and dotted lines indicate the de Vaucouleurs profile, projected Jaffe model, 
exponential profile, and our approximated exponential model, respectively, (b) Surface mass within R. Lines are the same as in (a). Note that our approximated 
exponential model agrees fairly well with the exponential profile. 



p{r) 



dfi{R) dR 



(A3) 



However, since the obtained profile is generally not given by an analytic function, it is too complicated for analyses to use the exact profile 
Therefore we adopt a simple, analytic model approximating the de Vaucouleurs profile when projected, that is, the Jaffe model (laffe 1983) 
By introducing characteristic density pt and size ri, of baryonic matter, the density distribution p(r) is written as follows, 



p(r) = 



(A4) 



where the half-mass radius in three-dimensional distribution coincides with r^. The normalization of the density distribution is given by 
equating the total masses of the de Vaucouleurs and the Jaffe models. From this condition, we find 



2520 exp(srdv)-Re 

Pb = g 3 P-e 



(A5) 



Moreover, we also find that the relationship between three- and two-dimensional half-mass radii, and Re, is Re — 0.744re. These relations 
lead to ^e ~ 4pi,re. 

The profile of surface density p{R) and the profile of surface mass AhD{< R) within R are shown in Figure Al. The solid and dashed 
lines indicate the de Vaucouleurs profile and projected Jaffe model, respectively. We adopt the normalizations of density and size given by 
the above equations. As proved by many authors, the Jaffe model well reproduces the de Vaucouleurs profile. 



A2 Approximate density distribution for the exponential profile 



Spiral galaxies and dwarf elliptical galaxies have an exponential luminosity profile, 
R 



p{R) = pe exp 



(?cxp 



(A6) 



where poxp = 1.68. Assuming the spherical symmetry, the deprojected three-dimensional density distribution using equation (A3) is ex- 
pressed by the modified Bessel function of the second kind, ^effoxp exp{gcxp)Ko{gc^pr / Re) /irRe- However, because of simplicity, we use 
an approximate form for three-dimensional density distribution. Note that the modified Bessel function of the second kind Ko is asymp- 
totically proportional to exp{—gcxpr / Re)/ \/r at r oo. Then, to obtain a simple form expressed in elementary functions, we adopt the 
following density distribution 
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p{r) = Pb\ —exp 



^-1 
rt 



As mentioned in the previous subsection, equating the total masses of these distributions, we find 



Pb 



rPe, 



(A7) 



(A8) 



and _Re — 0.978ri,. The half-mass radius for three-dimensional distribution is given by re — 1.295ri,, then Re — 0.755re. The projected 
density and mass profiles are also shown in Figure kl The dot-dashed and dotted lines indicate the exponential profile and our approximated 
exponential model, respectively. It is found that the exponential profile is well reproduced by our model. 



APPENDIX B: DERIVATIONS OF DYNAMICAL RESPONSE IN SPECIFIC CASES 

In this section we derive the relation between density and size for the adiabatic change. [The instantaneous change can directly be solved by 
substituting the function f{z) into equation M.] 



Bl de Vaucouleurs-like distribution of baryon in the isothermal dark halo 

Using the density and mass distributions of baryon and dark matter, the self- and interaction-potential energies of baryon in the isothermal 
halo of dark matter are calculated by 



USn^Gplrl 



-WavAso = G4-K-^Gpbpdrlrl 



ln(l + z) 



In ( 1 + - 

2 



respectively. From equation (g), the differential equation for adiabatic gas removal is 

dy ^ _ ln(l + z) 

dz z 2z^ ' 



then we obtain a solution 



y = 



c 



1 

2z4 



[-z + {l + z)ln{l + z)], 



where C is an integration constant given by equation (|l3[). In terms of mass instead of density, the response can be expressed by 

m = 4yz^, 

where m denotes the mass ratio of baryon to dark matter. 



(Bl) 
(B2) 

(B3) 
(B4) 
(B5) 



B2 de Vaucouleurs-like distribution of baryon in the homogeneous dark halo 

The interaction-potential energies of the de Vaucouleurs-like profile of baryon in the homogeneous halo of dark matter is 



TTrint 647r 3 2 

H/dv,hom = -^;—Gpbpdrbrd 



The differential equation for adiabatic gas removal is 
1 



dz z 



6z'^ 



2_i + Mi±£) 

z z^ 



2zln 



hi) 



Then we obtain 

- £. _L 

zi'^ 24^4 

where C is an integration constant. In terms of mass, we obtain 
m = 12j/2^. 



62 + - 22^ + 22* In (^1 + i ) - 2(3 + 22) ln(l + 2) 



(B6) 



(B7) 



(B8) 



(B9) 



B3 Exponential-like distribution of baryon in the isothermal dark halo 

The self-potential energy for baryon is given by 



■k^{3tv - 4). 



=29 



-GpbTb , 



(BIO) 
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and the interaction-potential energy for baryon embedded in the dark halo is 



-W^eSp.iso = Gpbpdrlrl- 



2gexp{-g/z) 



+ {37 + ^ - (5 + ^) erf (y^) + ln64 + 31n g) } + SV^G^ ( 



29 



From this, the differential equation becomes 



dy y 



2jZL£+6^exp(-f)+V^< 



5 5 

then integrating the differential equation, we obtain 



3/2 ' 



erf 



7/2„-9 



exp 



(-!) 



3/2 



(37r-4)^;4 

where C is an integration constant. In terms of mass, we obtain 

30FeS 3 



B4 Exponential-like distribution of baryon in the homogeneous dark halo 

The interaction-potential energy for baryon embedded in the dark halo is 



T,rint ^ 3 2^6" 

l^exp.hom = (^PbpdTbra 



Therefore the coefficient b is 



--(-f){«v^-™y!--(0")-'^( 



105% -36-1-16- 



)erf( 



Qg5/2 

and the function f{z) is 

«^)-^— (-1) 

The differential equation for y is 



^ - 70, / - - 105 - 
^ V 9 \9 



3/2 



V^(l05£;-36 + 16f)erf(^) 



dz z 3{3n - 4:)z'^ 



exp I -M j 56 140^ +21qI^ 



8^/n Vtt < 8 - 105 ( - 1 ^ erf 



then integrating the differential equation, we obtain 



C 



^ z4 ^ 6(37r - 4)z4 



exp 



(-1) 



-24 + 28- +70^ -M05^ ] -|-8V7r(3-2^ 



erf 



(7!) 



48 -32- -I- 105^ 



where C is an integration constant. In terms of mass, we obtain 

9y^eO 3 



4g5/2 



■yz 



B5 de Vaucouleurs-like distribution of baryon in the NF W dark halo 

The NFW density distribution is written as 



p(r) = PdC^ 



where c is the concentration parameter. The mass within radius r is 



V TdJ 1 



cr/rd 



+ cr/rd 
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(B22) 



The interaction-potential energy is 



dV,NFW 



eATvGpbpdTtrd 



cln(2) ln(l — cz) In(ca) + Li2(cz) 



1 — cz 



where Li2(z) is the dilogarithm, J°[ln(l — t)/t]dt. The differential equation is 



^ = + — 

dz z 2z^ 



1 — cz (1 — C0)2 
then integrating this equation, we obtain 



cfl — 2c2:) , . , ln(l — cz) ln(c2:) lA^icz^ 
m[cz) — — ^ — - 



y 



C 1 



2cz4 



-Acz + (4 - cz) ln(l - cz) ln(c«) + """"^^ ^ 3cz)^ln(cz) _^ _ ^^)Li2(cz) 



1 — cz 



where C7 is an integration constant. In terms of mass, we obtain 

-,3 



m = 



ln(l + c) -c/(l-|-c)' 



(B23) 



(B24) 



(B25) 



(B26) 
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